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The first test of the Kugo-Ojima colour confinement criterion by the lattice Landau 
>• I gauge QCD simulation is performed. The parameter u which is expected to be —16^ in 
the continuum theory was found to be —0.76^ in the strong coupling region. The data is 
O ■ analysed in connection with the theory of Zwanziger. In the weak coupling region, the 
^ . expectation value of the horizon function is negative or consistent to 0. 

O _ 

^ . 1. Introduction 

There are various manifestation of colour confinement in QCD. One is the linear poten- 
Ph! tial between quarks, which appears in the quenched lattice QCD simulation. A mechanism 
^ I for the appearance of the linear potential was proposed as well by Gribov [|^ about 20 
years ago. He showed that the Landau gauge fixing or the Coulomb gauge fixing does not 



fix the gauge field uniquely and that the restriction of to a physical space which 
is called the Gribov region enhances the singularity in the ghost propagator and induces 
^ ', the linear potential. 

In the field theory, the confinement implies absence of free single coloured particle state 
in the asymptotic Hilbert space. The physical Hilbert space should satisfy symmetry 
specified by the Lagrangian, and the QCD Lagrangian is invariant under the BRS(Becchi- 
Rouet-Stora) transformation. Kugo and Ojima[0 proposed in 1978, a criterion for the 
colour confinement in the Landau gauge, based on the BRS symmetry, which consists of 
a two-point function produced by the ghost, the antighost and the gauge field becomes 
—S^, where a and b specify the colour in the adjoint representation. Analytical calculation 
of this value is extremely difficult and so far no verification was performed. 

In the lattice QCD, the Gribov region still does not define the gauge field uniquely but 
there is a unique minimum in the fundamental modular region]^. He argued that the 
restriction to the fundamental modular region implies the regularity of the horizon tensor, 
whose transverse projection is identical to the — of the two-point function of Kugo-Ojima. 

Both the Gribov-Zwanziger's theory and the Kugo-Ojima's theory suggest that in the 
lattice Landau gauge, the gluon propagator is infrared finite, which is confirmed by the 
lattice QCD simulation [H|-B]. In this paper we simulate the two-point function specified 
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by the Kugo-Ojima in the lattice Landau gauge. 

2. The Kugo-Ojima confinement criterion and the Gribov-Zwanziger's theory 

2.1. Kugo-Ojima's theory 

A sufficient condition of the colour confinement given by Kugo and Ojima[0] is that 
defined by the two-point function of the FP (Faddeev- Popov) ghost fields, c(x), c(?/), and 

My), 

I e^Pi-y)l^Q\TD,c\x)g{A, x c),(i/)|0)rfx = {g,, - ^K{p') (1) 
J p 

satisfies ^^(O) = ~5^. 

Essential points in their argument is based on the BRS invariance of the QCD La- 
grangian accompanied by the gauge fixing and FP terms. 1) The Ward-Takahashi iden- 
tities implies that the gauge fields the auxiliary field 5", the covariant derivative 
of the ghost field D^c"' and the antighost field c"(x) necessarily have massless asymptotic 
fields which forms the BRS-quartet. The Hilbert space is decomposed into the FP ghost 
number eigenstates and in this bases, the BRS-quartet space have the zero-norm. 2) The 
BRS charge Qb is conserved, and the Noether current corresponding to the conservation 
of the colour symmetry is 

9j; = d^F;, + {QB,D,c}, (2) 

where its ambiguity by divergence of antisjTiimetric tensor should be understood, and 
this ambiguity is utilised so that massless contribution may be eliminated for the charge, 
Q"", to be well defined. The massless component in the current {Qb,D^c} is absent if 
1 + u = 0. 

2.2. Gribov-Zwanziger's theory 

The Landau gauge of the QCD is specified by (9^v4^ = 0. The Gribov region Q is 
specified by ensemble of local minimum points of gauge orbits under the variation with 
respect to (? = as follows. 

A\\A^f = -2{dA\e) + (e| - dV\e) + ■■■ (3) 
Q = {A\-dV>0 , dA = 0} . (4) 

The physical space of the gauge field is characterized by the condition that the FP 
determinant is positive. In the Coulomb gauge, the singulatity of the ghost propagator 
yields enhancement of the infrared singularity of the Coulomb potential 

In the lattice simulation, the unique gauge field configuration can be attained by the 
restriction to the fundamental modular region A^], which is specified by the absolute 
minimum along the gauge orbits. 

A = {A|pf = MiUgllA^f }, Acq. (5) 

Let the gauge configuration be in the fundamental modular region obtained by the 
optimising function I{U^) = y^(l Retrll^ ), and let it be a global minimum even 

X,IJ. 
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under the gauge transformation of larger period (the region is called the core region), and 
the two point tensor be defined as G^^xy^'^^ = (tr ^A"^-D^ — 'dD^~^'^^'^^^ 
Zwanziger's theory[^,10], the horizon function H{U) defined as 



xy 



liniG^^(p) - e 



(6) 



is negative or in finite volume, and in the infinite volume limit. Here N is the number 

of colour, V is the lattice volume, and e = —RetrUi) /V . 

I ^ 

Note that in a, d dimensional lattice e/d = 1 if all links = and that the value of 
1 — e/d has the meaning of the distance from this vacuum. 

3. The Lattice simulation 

We define the gauge fieldf^ on links as an element of SU{3) Lie algebra as, and perform 
the gauge transformation as e^^'^^ = gle'^'^-'^gx+ij. 

The Landau gauge is realised by minimising ||^^||^ via a gauge transformation g'^Ug, 
where = e^. e is obtained by solving the equation with a suitable parameter rj 

Me = -d^D^{A)e = 7]dA (7) 

The obtained norm ||y4|| is close to that obtained after the smeared gauge fixing|]^ 
within 1%. 

The inverse FP operator, 7Vl~-'^[f/] = (Mq — Mi[U])~^, is calculated perturbatively by 
using the Green function of the Poisson equation[Q. 

In use of colour source |A°a;) normalised as tr{\"'x\X'^Xo) = S"'''6x,xoy the ghost propagator 
is given by 

C^'ix^y) = {tr{X'^x\iM[U])-'\X'y)) (8) 

where the outmost () specifies average over samples U. 

The ghost propagator is infrared divergent and its singularity can be parametrised as 
Z3p-2-a = Z3(4sin2 f The data of /? = 5.5, 8^ x 16 lattice indicates that the 

singularity is approximately while the data of 16'^ lattice indicates that it is p~^-^. 

The finite-size effect is not so large j^. These qualitative features are in agreement with 
the analysis of the Dyson- Schwinger equation[0. 

Using the above inverse FP operator, we obtained that ^^(0) is consistent to —c6^, c = 
0.7 in SU{3) quenched simulation, /? = 5.5, on 8"^, 12^ and 16"^. Fig. 2 shows the value of 

In terms of the Kugo-Ojima parameter c, the function 6*^^(0) can be written as 1 + 3c. 
In the 16"^ lattice, e/d = 0.78 and 0.86 for /? = 5.5 and P = 6, respectively. They are 
numerically close to the inverse of the renormalization factor of the ghost propagator 
I/Z3. The Slavnov- Taylor relation Zi = Z3/Z3 = 1 implies e/c? ~ I/Z3. 

Our data of 16*^ lattice suggest that when P becomes larger, 6*^^(0) becomes smaller, 
while the Zwanziger's parameter e has the opposite tendency. 

It is to be remarked that in the Zwanziger's theory, Kugo-Ojima criterion c = 1 does 
not hold in view of e/d 7^ 1. In the case of /5 = 5.5, our data of Gfj^fj^{0) agrees with e. 
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Figure 1. The ghost propagator as func- 
tion of the lattice momentum. The data 

are /? = 5.5, 16"^. The fitted curve is 
1.287/p2-779. 



Figure 2. The Kugo-Ojima parameter \u'^\ 
as the function of the spatial extent of 
the lattice aL{fm). The data are j3 = 
6, 8^ 12^ 16^ and (3 = 5.5, 8^ 12^ 16^ from 
left to right, respectively. 



Jl c G^^jO) e e/d S I/Z3 

5.5 0.712(18) 3.14(5) 3.13(1) 0.783 0.657 0.777 

6.0 0.628(56) 2.88(17) 3.45(1) 0.863 0.694 0.860 



However, when the optimising function I{A^) = = '^trA^^A^ is used instead of 



X,fJ, 



I{U^), the function e/d is to be replaced by the link average {S{A)) = ( ^^1^(4/2 )' where 
A = adjA, and the value is reduced by about 20%. The positive horizon function for 
f] = 5.5 implies our configurations are not in the core region. In the weak coupling region 
{P = 6.0), the horizon function is negative or consistent to 0. 
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